arXiv:1504.05060vl [cond-mat.str-el] 20 Apr 2015 


Symmetry of the Fermi surface and evolution of the electronic structure across the 
paramagnetic-helimagnetic transition in MnSi/Si(lll) 

Alessandro Nicolaou,^’* Matteo Gatti,^’^’^ Elena Magnano,^ Patrick Le Fevre,^ Federica Bondino,^ 
Francois Bertran,^ Antonio Tejeda,^’^ Michele Sauvage-Simkin,^ Alina Vlad,^ Yves Garreau,^’^ 
Alessandro Goati,^ Nicolas Guerin/ Fulvio Parmigiani/’^’and Amina Taleb-Ibrahimi^ 

^Synchrotron SOLEIL, L’Orme des Merisiers, Saint-Aubin, BP 4S, F-91192 Gif-sur-Yvette Cedex, France 
^Laboratoire des Solides Irradies, Ecole Poly technique, CNRS-CEA/DSM, F-91128 Palaiseau, France 
^European Theoretical Spectroscopy Facility (ETSF) 

^ CNR-IOM, Laboratorio TASC, S.S. I 4 Km 163.5, 3414^ Basovizza, Trieste, Raly 
^ Universite de Lorraine, UMR CNRS 7198, Institut Jean Lamour, BP 239, F-54506 Vandoeuvre-ls-Nancy, France 
^MPQ, Universite Denis Diderot Paris VII, Bat. Condorcet, F-75205 Paris Cedex 13, France 
'^Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, 3012 Bern, Switzerland 
^Elettra-Sincrotrone Trieste, Strada Statale I 4 , 3414^ Basovizza, Trieste, Italy 
^ Universitd degli Studi di Trieste, Via A. Valerio 2, 34127 Trieste, Italy 
^^International Faculty - University of Cologne, Germany 
(Dated: August 16, 2016) 

MnSi has been extensively studied for five decades, nonetheless detailed information on the Fermi 
surface (FS) symmetry is still lacking. This missed information prevented from a comprehensive 
understanding the nature of the magnetic interaction in this material. Here, by performing angle- 
resolved photoemission spectroscopy on high-quality MnSi films epitaxially grown on Si(lll), we 
unveil the FS symmetry and the evolution of the electronic structure across the paramagnetic- 
helimagnetic transition at Tc ~ 40 K, along with the appearance of sharp quasiparticle emission 
below Tc. The shape of the resulting FS is found to fulfill robust nesting effects. These effects 
can be at the origin of strong magnetic fluctuations not accounted for by state-of-art quasiparticle 
self-consistent GW approximation. From this perspective, the unforeseen quasiparticle damping 
detected in the paramagnetic phase and relaxing only below Tc, along with the persistence of the 
d-bands splitting well above Tc, at odds with a simple Stoner model for itinerant magnetism, open 
the search for exotic magnetic interactions favored by FS nesting and affecting the quasiparticles 
lifetime. 

PACS numbers: 79.60.-i,71.20.Lp,75.50.Cc,71.27.+a 

Transition-metal monosilicides, such as TM-Si (TM = 

Mn, Fe, Go) and their solid solutions, show an intrigu¬ 
ing evolution of the electronic and magnetic properties 
as a function of temperature and doping [1, 2]. Among 
these materials, MnSi is a paradigmatic case for the ex¬ 
otic magnetic properties not accounted for by standard 
models of magnetism. At Tc 40 K in thin films) MnSi 
shows a phase transition from a paramagnetic metal state 
to a helimagnetic order. The effective magnetic moment 
yUe// of 2.27/iB/Mn drops to fisat = Q.d/r^/Mn in the 
saturated ferromagnetic phase [3, 4]. This value is sig¬ 
nificantly smaller than the l/r^/Mn value predicted by 
local-density approximation (LDA) calculations [5]. Fur¬ 
thermore, the magnetic transition can be progressively 
suppressed by applying a relatively low hydrostatic pres¬ 
sure of 14.6 kbar, while Tc tends to zero and an abrupt 
change of the resistivity from the T^ Fermi liquid be¬ 
havior to the T^/^ non-Fermi liquid (NFL) character is 
observed. This NFL phase is quite robust with respect 
to temperature, pressure and magnetic fields, and it is 
accompanied by a partial magnetic order for p<21 kbar 
[6-8]. These properties are only partially accounted for 
by the Moriya’s self-consistent renormalisation theory for 
spin fluctuations [9], hence the nature of the magnetic 


transition and the mechanism of the strong correlations 
in MnSi is matter of contests [10, 11]. In this complex 
and fascinating scenario MnSi has been extensively stud¬ 
ied for the last five decades. Nonetheless, the unknown 
valence-band structure and FS symmetry have banned 
to unveil the nature of the magnetic interactions in this 
compound. The nesting of the Fermi surface (FS), as pro¬ 
posed by O. Narikiyo [12], could account for some of the 
anomalies observed. This topological property of the FS 
could indeed create correlation effects accompanied by an 
important damping of the quasiparticles (QP). Moreover, 
the strong FS instabilities produced by a perfect nesting 
could induce low temperature phase transitions in which 
ordinary Fermi liquid behavior is recovered [13]. Hence, 
the detailed knowledge of the FS topology is essential 
for a deeper understanding of the interplay between the 
complex magnetic phase diagram of MnSi and the NFL 
phases. Furthermore, a similar mechanism, in connection 
with the chiral character of the magnetic interactions, is 
expected to govern the properties of skyrmion lattices, 
discovered in this compound a few years ago [14]. Un¬ 
fortunately, the lacking of a cleaving plane in MnSi sin¬ 
gle crystals has prevented, so far, momentum-resolved 
photoelectron spectroscopies, prompting ns to epitaxially 



2 


grow very high quality MnSi layers on Si(lll). Details 
on sample preparation and data acquisition are reported 
in [15]. 

By performing angle-resolved photoemission spec¬ 
troscopy (ARPES) experiments on 5.5 nm thick MnSi 
epitaxial films, we have been able to disclose the FS sym¬ 
metry, revealing strong deviations of the experimental 
data from state-of-art band structure calculations as well 
as the potential for strong nesting for the most of the 
FS sheets, accompanied by significant QP broadening. 
Moreover, following the evolution of the electronic struc¬ 
ture across the paramagnetic-helimagnetic phase transi¬ 
tion at Tc ^ 40 K we detect, at some particular of 
the nested band, a notable sharpening of the quasiparti¬ 
cles. From this perspective, the unusual QP broadening 
of the paramagnetic phase which is relaxing only below 
Tc, along with the persistence of the d-bands splitting in 
the paramagnetic phase, in contrast with a simple Stoner 
model for itinerant magnetism [16, 17], point to novel 
magnetic interactions not accounted for by state-of-art 
quasiparticle self-consistent GW approximation (QSGW) 
[18, 19]. Hence, our data are disclosing a possible major 
role of the FS nesting properties for developing strong 
magnetic fluctuations affecting the quasiparticles lifetime 
and thus expected to govern the magnetic interactions in 
this compound. 

Figures l(a-b) show FS mappings of the FMR plane 
in the 2.17 to 5.16 kz range, across the second and 
third Brillouin zones (BZ), centered at r2=(0,0,2.38) and 
r3=(0,0,4.76) [20]. Both measurements were per¬ 

formed in the paramagnetic phase on two different sam¬ 
ples of the same thickness, but with different photon po¬ 
larization (linear horizontal LH and vertical LV). The FS 
in (a) was measured at 300 K in LH. The FS in (b) was 
obtained in LV at lower temperature (55 K) for reducing 
the thermal broadening of the quasiparticles. LDA and 
QSGW calculations results give the same description of 
the five 3d bands crossing allowing to use the com¬ 
putationally cheaper LDA for the FS analysis, however 
QSGW noticeably improve the agreement with the ex¬ 
perimental data for the Mn 3d - Si sp hybridized states 
at higher binding energies [see Fig. 2(f)]. In the bottom 
half of Fig. 1 (a) we plot the original LDA contours. In 
the top half of Fig. 1 (a), and in (b), we present the 
best match to our data, which was obtained by displac¬ 
ing each calculated band (B 1 to B5) with respect to 
by the amount indicated in the top row of Fig. 1 (c). We 
note that only by operating such a charge redistribution 
within the 3d manifold we can properly describe the ex¬ 
perimental data. Actually, while the LDA predicts three 
hole-like pockets centered at T, and two electron-like ones 
at R, we find instead that the biggest hole-pocket B3 is 
forming long tubes containing the R points and that the 
biggest electron-pocket B4 has increased in size. Con- 
trarily to B3, B4 is clearly visible also in the third BZ in 
LV polarization, where its FS contour and the one of Bl, 



FIG. 1: (a-b): Experimental FS mappings of the FMR plane 

done in normal emission and calculated using an inner poten¬ 
tial Vo = 7 eV [20]. (a) Intensity map at Ef of the second 
derivative of the EDMs from a 15-45 eV energy scan in LH. 
(b) Intensity map at Ef from a 50-98 eV energy scan in LV. 
The Ff values, plotted as gray circles in (b), were obtained 
by the Lorentzian fit of the MDC lines at Ef- Black arrows 
are the nesting vectors, (c) FS as predicted by the LDA with 
the original position of Ef (bottom) and with Ef shifted for 
the different bands in order to fit experimental data (top). 
Robust values in white, less compelling ones in gray, see text. 


recognized by its rhomboidal shape, were tracked more 
precisely by a Lorentzian fit of the momentum distribu¬ 
tion curves (MDGs) [gray circles in Fig. 1(b)]. This latter 
measurement reveals indeed that the hole-like pocket Bl, 
centered at F 3 and not visible in the second BZ and/or 
with LH polarization, has also increased in size. The 
certain identification of the B1,B3,B4 bands is provided 
by their characteristic shape together with the very good 
agreement with the intensity maps of Fig. l(a-b), mak¬ 
ing the corresponding energy shifts of 105, 230 and -155 
meV, quite robust values. Less compelling are the values 
proposed for B2, lifted up by 390 meV to better account 
for the intensity parallel to the MR directions, and for 
B5, not shifted, as the present data do not allow to un¬ 
ambiguously identify them. 

In a similar way, as cited by Carbone in [ 21 ], S. B. 
Dougale and J. Laverock imposed a common rigid shift 
of 200 meV to the calculated bands in order to describe 
their positron annihilation measurements. However, this 
approach does not conserve the number of particles. Con¬ 
versely, to describe our data we have lowered Ef only for 
the three LDA bands B1,B2,B3 while we have increased 




















3 



-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 -1.0 -0.5 0.0 0.5 1.0 

kx(A'^) Binding en. (eV) k,(A^) 



3.5 4.0 4.5 5.0 5.5 -1.0 -0.5 0.0 M T R 

) Binding en. (eV) 


FIG. 2: (a,d) EDM and (b,e) corresponding EDCs of the 

FM and FR directions, acquired at T=55 K using LV polar¬ 
ization. The superimposed QSGW bands are shifted in energy 
to match the experimental FS of Fig. l(a,b). (c) Projection 
of the LDA bands of the FMR plane along FM. (f) Band- 
structure calculations in LDA (red solid lines) and QSGW 
(blue markers; linearly interpolated dotted lines are a guide 
for the eyes). 


it for B4, in a way that might conserve the total num¬ 
ber of particles. Even if a precise determination of the 
electron counting is meaningless without a complete 3D 
FS mapping, the observed disagreement with theory is 
unlikely originated from an uncontrolled doping of the 
sample and/or surface polarity effects, that would pro¬ 
duce the same shift of E^? for all bands, whereas it calls 
for other explanations. 

Interestingly, following this unforeseen charge redistri¬ 
bution, the experimental FS results more prone to nest¬ 
ing phenomena. The B1 FS sheet has indeed acquired 
a rhomboidal shape with parallel sections connected by 
nesting vectors 0.53 long, lying parallel to the [201] 
and [021] directions. Moreover, the top (bottom) half 
of B1 and the B4 band at lower (higher) k^ are almost 
perfectly nested, creating a coupling between hole- and 
electron-pockets with wave vectors 0.73 A“^ long. It is 
remarkable that this latter lies parallel to the [111] di¬ 
rection, along which the direction of the helimagnetic or¬ 
der is pinned below T^, and that, starting at about the 
same wave vector, neutron scattering experiments detect 
a strong temperature dependence of the magnetic exci¬ 
tations, not predicted by the Moriya’s spin-fluctuations 
theory [21]. Moreover, FS nesting is expected to produce 
unusual variations of the optical reflectivity through the 
enhancement of QP scattering [22], being thus a plausible 
responsible for the non-Drude optical behavior of MnSi 
reported in [11]. 

The indication that we are in presence of a correlated 


metal arises from the strong damping of the quasipar¬ 
ticles observed along both the FM and the FR high- 
symmetry directions, as presented in Fig. 2(a,b) and 
(d,e) respectively. On the top of the energy dispersion 
maps (EDMs) we plot the result of QSGW calculations, 
where the Mn 3d bands are shifted in energy according 
to the values used to reproduce our experimental FS. If 
the Fermi level crossings are clearly discernible from the 
FS maps, below Ep the lineshape becomes very broad 
and, accordingly, a significant mass renormalization is 
clearly visible in the top panel of Fig. 3(a), where the 
dotted line represents the B1 LDA band renormalized by 
a factor 3.3, a value which is consistent with previous 
Haas van-Alphen and infrared spectroscopy experiments 
[11, 23]. We note that while the LDA eigenvalues in 
principle cannot be interpreted as the electron addition 
and removal energies which are measured by ARPES, the 
renormalization and the rigid shift of the d-bands binding 
energies with respect to those calculated by the many- 
body GW approximation (GWA) [24], where the self¬ 
energy is given by a product of the one-particle Green’s 
function G and the screened Coulomb interaction W and 
therefore contains electron-hole and plasmon excitations, 
points to exotic magnetic interactions. Actually, the re¬ 
sulting dynamical screening of the photoemission hole is 
a crucial ingredient if we want to compare ARPES with 
band structure calculations, however magnetic correla¬ 
tion effects due to the coupling with spin excitations are 
not accounted for by the GWA [25-29]. Therefore, we 
conclude that this large broadening is intrinsic to the 
electronic excitations of MnSi and it is the fingerprint of 
a many-body mechanism associated with the spin degree 
of freedom. 

Even if part of this broadening arises from the super¬ 
position of the FS sheets corresponding to MnSi(lll) 
left- and right- handed twinned domains revealed by our 
GIXD measurements [15], this effect is expected to be 
smaller with respect to the observed one. Moreover, sev¬ 
eral factors, like final-state effects for 3D electronic struc¬ 
tures [30] or corrections to the sudden approximation 
[31], can make the measured photocurrent significantly 
differ from the intrinsic spectral function of a material. 
Nevertheless, those effects are expected to be strongly 
photon-energy dependent, while we do not observe any 
significant variation in 15-100 eV energy range. In ad¬ 
dition, we detect another spectral feature centered at F 
with the top at about 500 meV below Ep and not pre¬ 
dicted by GW spectral-function calculations [15]. Even 
though the dispersion of this feature in the FR direction 
is modest, it appears in a position forbidden for pure 
surface states, as it is already occupied by bulk bands, 
projected in Fig. 2(c) on the FM direction. Moreover, 
it does not fit the periodicity of the FM surface direc¬ 
tion, which is half of FM periodicity. The origin of this 
spectral feature, either a real hole-like band sunk below 
Ef or the incoherent component associated with one of 
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FIG. 3: (a) EDM images as a function of temperature ac¬ 

quired in LV with photons of 75 eV. They correspond to a cut 
of the reciprocal space parallel to FM but at kz = 4.5 
see dotted line in Fig. 1(b). (b) Comparison of the highest 
(dotted) and lowest (thick) temperatures EDCs of the data 
presented in (a), (c) EDCs as a function of temperature at 
ka; = 0.25 A“^. (d) Temperature dependence of EDM images 
taken at the same kz as (a) but rotated by 30 ° with respect 
to FM. (e) Fraction of the spectral intensity at the sides {kx > 
0.38 A“^) and between -0.1 eV and Ef as a function of tem¬ 
perature. In (a) and (c) solid lines represents LDA bands, 
and the dotted one LDA bands renormalized by 3.3. 


the Mn 3d bands at Ei?, calls for a deeper investigation 
which is beyond the scope of the present work. 

The effect of the magnetic transition on the electronic 
structure B1 is presented in Fig. 3. In (a) we show a cut 
of the reciprocal space parallel to the FM direction but at 
a lower [see dotted line in Fig. 1(b)]. Below 40 K we 
observe the appearance of sharp quasiparticle emission 
at about ka^=0.25 A“^, see Fig. 3(c). This sharp feature 
proves that the intrinsic limitation of the quasiparticle 
lifetime in the paramagnetic phase is relaxed when the 
system, cooled down below Tc, enters the ordered heli- 
magnetic phase. In Fig. 3(b) we compare the energy dis¬ 
persion curves (EDCs) for the highest and lowest temper¬ 
atures measured here, to show that the rising of the sharp 
feature below T^ unambiguously originates from an in¬ 
creased quasiparticle lifetime and not by a mere increase 
of energy resolution at low T. To complete the interpreta¬ 
tion of the magnetic transition we have precisely tracked 
the evolution of the electronic structure across T^. To 
this end, temperature-dependent cuts of the reciprocal 
space at the same k^^ but rotated 30° with respect to 
the FM direction, have been considered. Fig. 3(d) re¬ 
ports such data. By decreasing the temperature we de¬ 


tect the appearance of two symmetric bands at about 
±0.75 A“^. The evolution of these features through the 
magnetic transition has been tracked by following their 
relative spectral weight with respect to the whole EDM 
in Fig. 3(e). They have a constant weight below Tc, 
while at Tc they start to loose progressively weight until 
75 K, which is the maximum temperature measured here. 
It is known that the magnetic transition is suppressed in 
thin films by reducing the film thickness [32, 33]. The 
observed Tc and the thickness of our film measured by 
SQUID and GIXD respectively [15], are consistent with 
these studies, thus endorsing our novel recipe for grow¬ 
ing MnSi on Si(lll). Despite the fact that the spectral 
weight evolution of these bands allows to track correctly 
the temperature dependence of the magnetic transitions, 
they do not correspond to the calculated electronic bands 
and their identification calls for further investigations be¬ 
yond the scope of the present work. Nevertheless, the 
persistence of an effect of magnetism well above T^ im¬ 
plies the presence of local magnetic order well above Tq 
[16, 17], 

In conclusion, by performing ARPES experiments on 
high-quality MnSi epitaxial layers, we have disclosed the 
FS symmetry along with robust nesting effects that can 
be at the origin of strong electron-electron and electron- 
hole scattering, not anticipated by the state of the art 
of band structure calculations. In this scenario, the 
unusual quasiparticles damping detected in the param¬ 
agnetic phase and relaxing only below T^, along with 
the persistence of the d-bands splitting in the paramag¬ 
netic phase, in contrast with standard models of itinerant 
magnetism, point to novel magnetic interactions not ac¬ 
counted for by quasiparticle self-consistent GW approx¬ 
imation. Interestingly, this scenario seems to gain sup¬ 
port by following the evolution of the electronic struc¬ 
ture across the paramagnetic-helimagnetic transition at 
Tc ^ 40 K, where we observe an unforeseen sharpening 
of the quasiparticles at some particular kp values of the 
nested band. Henceforth, we propose here that the nest¬ 
ing properties of the MnSi Fermi surface, by affecting the 
quasiparticles lifetime through the development of strong 
magnetic fluctuations, govern the magnetic interactions 
present in this compound. 
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I. DATA ACQUISITION 

ARPES and GIXD measurements were performed at 
the beamlines GASSIOPEE and SIXS respectively, at the 
synchrotron SOLEIL. ARPES data were acquired using a 
Scienta R4000 hemispherical analyzer with the entrance 
slit orthogonal with respect to the plane containing the 
incident beam and the normal to the sample surface. The 
energy resolution for ARPES was set to 25 (35) meV and 
10 meV for the FS mapping in linear vertical LV polar¬ 
ization (linear horizontal LH) and the temperature de¬ 
pendences respectively. GIXD measurements were done 
using an UHV z-axis diffractometer at 18.41 KeV. 


II. ARPES DATA ANALYSIS 

In the presented ARPES data, each energy dispersion 
map (EDM) was normalized so that each EDG has the 
same area. The different EDMs where normalized one 
with respect to the other to have a continuous evolution 
of the intensity along the PR direction. This normaliza¬ 
tion procedure was used in order to reduce the effect of 
the matrix elements relative to the different final states 
for different the photon energies, and to highlight the 
symmetry of the FS surface in Fig. 1(b) and the disper¬ 
sion of the electronic bands in Fig. 2 (d-e) of the main 
text. 


III. STATE OF THE ART OF ARPES DATA ON 
MNSI 

Grytzelius et al. report the electronic structure of a 
very thin (5-8 A) MnSi film on Si (111). They inter¬ 
pret their result as the electronic structure of the Si(lll) 


surface in presence of MnSi terminations^. Kura et al. re¬ 
port ARPES measurements of an MnSi(lOO) single crys¬ 
tal polished and annealed; they observe some band dis¬ 
persion in the ligand Mn 3d Si sp states at high binding 
energy, but they do not resolve the Mn 3d bands and 
they do not observe any clear dispersion near Even 

higher energy resolved measurements of J.Y. Son et al. 
reported^ about a scraped MnSi single crystal, do not 
show considerable variation of the band structure with 
temperature and look like the ones obtained on an oxi¬ 
dized sample. We stress the fact that one of the main 
limitation of the ARPES technique is its surface sensi¬ 
tivity (limited to the very first atomic layers in the UV 
regime) and the necessity of an atomically flat surface, as 
the measure of the k//, the component of the wavevector 
k parallel to the sample surface, is done by measuring the 
take-off angle with respect to the sample surface. There¬ 
fore, surface preparation of single crystals might change 
the stoichiometry of the first atomic layers, preventing a 
measure representative of the bulk and, likely, the mea¬ 
sure of the magnetic transition which is strongly doping 
dependent in this compound. This might be the reason 
why, so far, no evolution of the electronic structure as a 
function of temperature across Tc has been observed. 


IV. DETAILS ON SAMPLE SYNTHESIS AND 
CHARACTERIZATION 

The epitaxial MnSi films were prepared by simulta¬ 
neous co-evaporation from Knudsen cell crucibles of Mn 
and Si on a Si(lll)-7x7 substrate at room temperature, 
and subsequent annealing at 250 °C for several hours. 
With respect to our previous works^’^, the long anneal¬ 
ing time at a relatively low temperature of 250 °C was 
found to be the key to obtain homogeneous films showing 
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a well-defined band dispersion. The absence of spurious 
phases is proved by the grazing incidence X-ray diffrac¬ 
tion (GIXD) and the presence of the magnetic transition 
by SQUIDS 


-1h @ 250 °C 


1h@ 250 °C -4h40@ 250 °C 



-4 -3 -2 -1 0-4 -3 -2 -1 0-4 -3 -2 -1 0-4 -2 0 

Binding Energy (eV) Binding Energy (eV) 



FIG. 1: (a-b) ARPES data after Ih heating at 250 °C in 

two different position on the sample, (c-d) ARPES data after 
4h40min heating on an homogeneous sample and correspon¬ 
dent stacks, (e) Photoemission data obtained by integration 
of EDMs in (a-c) over the whole k-range. All these data 
where taken at room temperature using a 30 eV LH polarized 
beam, (f) 1 ROD at h(k)=4 compared to the fit. Data and fit 
where vertically shifted for clarity, (g) h(k),l map of 5.5 nm 
MnSi(lll)/Si(lll). GIXD measurements where performed 
at room temperature and ambient pressure on Si capped 
MnSi/Si(lll) sample. h,k,l am expressed in the MnSi(lll) 
basis. The bright spots not highlighted are the ones of the 
Si(lll) substrate. 


In figure l(a-c) we show how the heating time strongly 
influences the homogeneity and crystallinity of our films. 
Higashi et alJ studied the temperature dependence of 
the FEED spot intensity and found that the best heating 
temperature is 250 °C, corresponding to the maximum 
spot intensity, and therefore a well defined long-range or¬ 
der. In their work, they heated for 30 min 3 monolayer 
of Mn evaporated on the Si(lll) substrate. Looking at 
our ARPES data, we have found that in order to ob¬ 
tain homogeneous samples showing a well defined band- 
dispersion it is necessary to heat at this relatively low 
temperature for a long time. The data presented in Fig. 
l(a,b) where both taken after 1 hour of heating at 250 °G, 
but in different parts of the sample. As is clear from these 
two data sets, after only 1 h at 250 °C the sample is not 


homogeneous and presents weakly dispersing bands only 
in some regions. Additional heating of the same sample 
for about 4 hours, yielded the nice dispersing bands pre¬ 
sented in l(c,d), which were found to be homogeneous 
over the whole sample. The resulting dispersion of these 
bands is clear by looking at the corresponding Energy 
Dispersion Curves (EDC) presented in 1(d). If we look 
at the comparison of the integrated photoemission data 
presented in 1(e), obtained by integrating each energy 
dispersion map (EDM) over the whole k-range, we ob¬ 
serve that, despite all of them present clear weight at 
the Fermi level (E^?), only in the long heated sample the 
Mn 3d states (max @ -0.12 eV) are well separated with 
respect to the bonding Mn 3(i hybridized with the Si sp 
bands, which are found at -1.64 eV. 

The homogeneity of the long heated sample and the 
presence of the single B20 crystalline phase is confirmed 
by the GIXD measurements presented in Fig. l(f-g). The 
diffraction pattern obtained by GIXD shows a very good 
epitaxy of the film with the crystallographic relationship 
with the substrate being MnSi(lll)/Si(lll). From re¬ 
ciprocal space maps parallel to the surface plane (not 
shown), we have identified the surface cells of the sub¬ 
strate and of the film, and we have measured a mismatch 
between the film and the substrate of 2.2%, meaning that 
our film is not fully relaxed at its own lattice parameter, 
but 1.1% tensile strain is still present^. We recall that 
a fully relaxed MnSi(lll) film on Si(lll) present a 3.3% 
difference in the in-plane lattice parameter with respect 
to the substrate. From the the h(k),l map of the recip¬ 
rocal plane in the base of MnSi presented in Fig.l (g) we 
have determined the the presence of 2 twinned domains. 
The two different sets of diffraction spots, expected from 
each B20 structure, are presented as circles and squares 
for the two domains. The presence of twinned domains is 
due to the growing of a non-centrosymmetric structure as 
the B20 one on a centrosymmetric substrate (Si(lll))^. 
Karhu et al. reported domains of few jj? size^. From 
a fit of our hk scans with a gaussian line shape we can 
estimate a minimum size of individual twin domains in 
our sample to less than 0.01 p?. ARPES measurements 
were performed with a photon beam of about 100x100 
size; we therefore illuminate both domains and we expect 
the mere superposition of the two band structures. By 
fitting the rod at h=k=4 presented in 1 (g) with ROD 
we can estimate the presence of 45 % of one domain and 
55 % of the other. We can also estimate precisely the 
sample thickness to 55 A. The best fit was obtained by 
stacking 7 unit cells of MnSi and leaving free the c-strain 
of each unit cell, yielding a strain variation from -1.75 % 
at the interface to 0.5 % at the surface. By performing 
EDA calculation using a trigonally stressed B20 unit cell, 
both in and out the <111> plane, we have observed only 
negligible variations of the electronic structure. 

In Eig. 2 (a) and (c) we present Eermi surface (ES) 
mappings acquired at fixed photon energies of 35 eV and 
75 eV using linear vertical and linear horizontal polar¬ 
ization respectively. These polar scans are compared 
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FIG. 2: (In (a) and (c) we present Fermi surfaces mapping 
acquired by performing polar scans at photon energies of 35 
eV and 75 eV respectively. These are compared to LDA band 
structure calculations where the electronic bands were shifted 
in energy to reproduce the experimental Fermi surfaces pre¬ 
sented in Fig. 1 of the main text. These two measurements 
probe planes of the reciprocal space parallel to the FMM one 
but shifted along the 111 direction of -16 % and - 67 % of the 
FR direction. The meshes used for the LDA calculation are 
presented in (d) together with the simple cubic Brillouin zone 
of the B20 crystal structure and the FMM plane (green) and 
the FMR one (red). In (b) we show the Low energy electron 
diffraction pattern of the MnSi thick him. 


with LDA band-structure calculations performed using 
an hexagonal mesh of the same size of the FMM plane 
and parallel to it, but shifted along the 111 direction of 
-16 % and - 67 % of the FR direction length. The calcu¬ 
lated bands were shifted in energy by the same amount 
used to reproduce the experimental FS obtained by the 
energy dependencies and presented in the Fig. 1 of the 
main text. The good agreement between data and this 
simple model also away from the FMR plane further con- 
hrms the need to readjust the energy of the bands cal¬ 
culated by LDA in order to describe the experimental 
Fermi surface of this compound. The low energy diffrac¬ 
tion pattern (LEED) is presented in (b), to show the 
long range order of the surface and no presence of sur¬ 
face reconstructions. The small hexagon correspond to 
the surface termination of the MnSi film, while the big¬ 
ger one to the MnSi bulk. The sixfold symmetry of this 
latter is likely due to the presence of twinned domains. 

The magnetic transition of our 5.5 nm thick 
MnSi/Si(lll) sample was characterized by means of 
SQUID magnetometer (MPMS Quantum Design), as 


presented in Eig. 3. The MnSi thin film plane was set 
parallel to the applied magnetic field H so that the in¬ 
plane long moment is recorded. In a first strategy, in Eig. 
3 (a), after a zero field cooling, H = 0.05 T was set and 
the moment was then measured for increasing tempera¬ 
ture. In this way we observe a transition at T^ = 40 
K, after the subtraction of the Si(lll) background. In 
a second strategy, in Eig. 3 (b), the sample was cooled 
down to 5 K and saturated by applying the field H = 5 
T. After removing this field, the remnant moment was 
measured upon warming, yielding again a T^ of 40 K. 
Eor all curves, we observe a continuous and straight de¬ 
crease starting from 5 K down to T^, which is compat¬ 
ible with weak itinerant ferromagnetism Above T^, the 
background values are still positive even when the applied 
field was totally cancelled. Thus the sample may show an 
additional ferromagnetic component (even if weak) up to 
350 K, which is maximum temperature measured here, 
and attributed in^^^i surface ferromagnetism. 




Temperature (K) 

FIG. 3: In-plane magnetic moment under the held H = 0.05T 
for the thin him MnSi(lll). Zero held cooling and subtraction 
of the Si(lll) background moment 
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V. THEORETICAL FRAMEWORK AND 
COMPUTATIONAL DETAILS 


In the GW approximation (GWA)^^ the (non-local, dy¬ 
namical and complex) self-energy E is given by the con¬ 
volution in frequency space: 


I](ri,r2,w) = ^ J duj'e"^'^'G{ri,r2,uj+co')W{ri,r2,co') 

( 1 ) 

between the time-ordered one-particle Green’s function 
G: 


G(ri,r2,w) = ^ — 


4>i{Yl)(f>*{Y2) 


uj — £i+ ir]sigii{ei — jj) ’ 


( 2 ) 


(where fi is the Fermi energy) and the dynamically 
screened Coulomb interaction W: 

W(ri,r2,co) = y dr3e“^(ri,r3,a;)u(r3 - r2). (3) 

Here v is the (static) bare Coulomb interaction and 
is the inverse dielectric function that describes screening 
of V through electron-hole and collective plasmon exci¬ 
tations. In the GWA is calculated in the random- 
phase approximation (RPA). Within the quasiparticle 
self-consistent GW (QSGW) scheme^^’^"^ the quasipar¬ 
ticle wavefunctions (j)i and energies Si are obtained from 
the (static and Hermitian) potential: 

^ ^ E Pb'fe) + (<^,-|. (4) 

At self-consistency the quasiparticle energies £i deter¬ 
mine the band structure and the spectral function is cal¬ 
culated from 

( 5 ) 

We have used the QSGW implementation of Ref. 
15 in the Abinit code^^. We have adopted the 
MnSi experimental crystal structure from Ref. 17. 
We have employed Troullier-Martins norm-conserving 
pseudopotentials^^ with Mn 35 and 3p semicore states 
treated as valence electrons and 160 Hartree energy cut¬ 
off. We have used a 6x6x6 F-centered k-point grid and 
350 (500) bands for the calculation of the screening (self¬ 
energy). The plane-wave size of the dielectric function 
corresponds to 5.3 Hartree cutoff. Finally, the frequency 
integration in Eq. (1) has been performed with a contour- 
deformation technique using 30 frequencies on the real 
axis and 6 on the imaginary axis. 


VI. CALCULATED BAND STRUCTURE AND 
SPECTRAL FUNCTION 

Fig. 4 shows the band structure of the paramagnetic 
phase of MnSi calculated in LDA and QSGW. LDA re¬ 
sults are in agreement with those reported by Jeong 



FIG. 4: Band structure calculated in LDA (red lines) and 
within QSGW (blue markers; linearly interpolated dashed 
lines are a guide for the eye). 


and Pickett in Ref. 19. Since Kohn-Sham eigenval¬ 
ues in principle are not the electron removal energies 
measured by ARPES and in order to take into account 
non-local exchange and correlation effects, we have used 
a Green’s-function approach within the GWA. Due to 
the the presence of localized Mn 3d states, the accurate 
QSGW scheme has been employed (beyond standard per¬ 
turbative corrections to LDA band structure^^). 

In the comparison with the LDA band structure, 
QSGW noticeably differs at energies below ^ -3 eV, in¬ 
creasing the energy separation between Mn 3d^ located 
around the Eermi energy, and Si 3p at higher binding 
energies. Si 35 bands start below -7 eV (not shown). 
The same differences between LDA and QSGW appear 
for the conduction states, where Si 3p are present above 
^ 1 eV. On the contrary, for the Mn 3d bands crossing 
Ep LDA and QSGW are overlapping almost completely. 
This justifies the use of the LDA for plotting the Eermi 
surfaces of MnSi. 

We have also calculated corrections to the LDA due 
to spin-orbit coupling and vertex corrections beyond 
the GWA (in the GWF approximation with an LDA 
vertex^^) without finding qualitative changes. Also in the 
helimagnetic phase (approximated as ferromagnetic in 
the calculations) the QSGW magnetic moment decreases 
only by 0.04 /r^/Mn with respect to the LDA^^, confirm¬ 
ing that the theoretical overestimation of the magnetic 
moment should be fixed by taking into account spin ex¬ 
citations not included in the GWA. Similar differences 
in the band structure between the LDA and QSGW as 
those discussed here for the paramagnetic phase are seen 
also for the magnetically ordered phase. 

The top panel of Eig. 5 shows the GW spectral func¬ 
tion for the top-valence and bottom-conduction states at 
F for paramagnetic MnSi. The quasiparticle peaks have 
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Spectral function at k=r 



Energy (eV) 


Total spectral function 



all a finite width corresponding to a finite lifetime. The 
incoherent spectral weight appearing below Ep that is 
associated to the conduction states (i.e. having quasipar¬ 
ticle peaks above E^?) is too small to explain the features 
seen by ARPES at ^ 0.5 eV binding energy (see main 
text). Einally, the bottom panel of Eig. 5 reports the 
GW spectral function summed over the valence and con¬ 
duction states and integrated over the whole Brillouin 
zone (i.e. it should be compared with angle-integrated 
photoemission experiments). 


FIG. 5: (Top panel) Spectral function for the top-valence 
and bottom-conduction states at F. (Bottom panel) Total 
spectral function summed over bands and integrated over the 
whole Brillouin zone. The zero of the energy axis is the Fermi 
level. 
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